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Abstract. We introduce quasi-local (QL) scalar variables in spherically symmetric LTB models. If the QL scalars are denned 
as functionals, they become weighed averages that generalize the standard proper volume averages on space slices orthogonal 
to the 4- velocity. We examine the connection between QL functions and functionals and the "back-reaction" term i? in the 
context of Buchert's scalar averaging formalism. With the help of the QL scalars we provide rigorous proof that back-reaction 
is positive for (i) all LTB models with negative and asymptotically negative spatial curvature, and (ii) models with positive 
curvature decaying to zero asymptotically in the radial direction. We show by means of qualitative, but robust, arguments that 
generic LTB models exist, either with clump or void profiles, for which an "effective" acceleration associated with Buchert's 
formalism can mimic the effects of dark energy. 
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INTRODUCTION. 

Recent observations apparently reveal that the universe is spatially flat and is undergoing an accelerated expansion. To 
account for these observations, a large variety of theoretical and empiric models have molded a dominant theoretical 
paradigm: the "concordance" model, based on the assumption that cosmic dynamics appears to be dominated by an 
elusive source ("dark energy") that behaves as a cosmological constant or as a fluid with negative pressure (see 0] for 
a review). 

The concordance model also assumes that inhomogeneities play a minimal role in cosmic dynamics at scales over 
100 Mpc, and thus can be adequately dealt with in terms of linear perturbations on a FLRW background. This 
assumption, and the existence of dark energy, has been challenged from various angles [2]. In trying to account 
with supernovae observations, numerous articles (see [2j for a review, see also Q) show that cosmic acceleration 
can be reproduced simply by considering large scale inhomogeneities in photon trajectories within the so-called 
"homogeneity scale" (100-300 Mpc). From a theoretical point of view, it has been argued that inhomogeneity implies 
that observations from distant highredshift sources must be understood in terms of averaged quantities 13 4 5,6,7,81, 
which in homogeneous conditions would be trivially identical with local quantities. Thus, non-linear spatial gradients 
of the Hubble expansion scalar and quasi-local effective energies would have an important effect in the interpretation 
of these observations |8 1. 

The spherically symmetric LTB dust models (9l [T0j [TTJ have often been used to test the effects of inhomogeneity in 
cosmic observations, as well as the issue of back-reaction in the context of Buchert's spatial averaging lfT2l[T3l . In the 
present article we examine the dynamical equations of there models in terms of suitably defined quasi-local variables 
and of spatially averaged scalars lfl4l [151 . We show that the back-reaction term is the difference between squared 
fluctuations of the expansion Hubble scalar, which in turn are related to spatial gradients of its average and quasi-local 
equivalent. Necessary conditions for a positive "effective" accelerations, mimicking the effect of dark energy, follow 
from comparing these fluctuations, and are satisfied as long as spatial curvature is negative in a sufficiently large 
averaging domain (even if smaller domains contain bound structures). Since these conditions are compatible with the 
observed void dominated structure, it is highly likely that such "effective" acceleration could be observed. However, 
further steps in this direction require a more elaborate numerical study of these models (see 1 16 1). 



LTB DUST MODELS IN THE "FLUID FLOW" DESCRIPTION. 



Spherically symmetric inhomogeneous dust sources are usually described by the well known Lemaitre-Tolman-Bondi 
metric and energy-momentum tensor in a comoving frame l9l [T0lfm 

ds 2 = -c 2 df 2 +^dr 2 + fl 2 (d0 2 + sin 2 ed</> 2 ). (1) 

T" b = pu a u b , (2) 

where R = R(ct,r), R' = dR/dr, & = J^"(r), u" = 8q and p(t,r) is the rest matter-energy density. The field equations 
G ab = KT ab (with K = 8 nG / c 2) for Q and ^ reduce to 

R 2 = 2 f +^ 2 -K (3) 

K 

2M' = KpR 2 R\ (4) 

where M = M(r) and R = u c 'V a R. The sign of J^" 2 — 1 determines the zeroes of R and thus classifies LTB models in 
terms of the following kinematic classes: 

3^ = 1 , parabolic models 

& > 1 , hyperbolic models 

< & < 1 , "open" elliptic models 

-1 < & < 1, "closed" elliptic models 

(5) 

where by "open" or "closed" models we mean the cases where the are, respectively, topologically equivalent to 
R 3 and S 3 . Regularity conditions ifTOlfTTIl require > for all r in hyperbolic models, while in open elliptic models 
can be, either negative for all r > 0, or it can have a zero 3^'(y) = 0, so that < for < x < y and > for 

x> y. 

Besides p and R given above, other covariant objects of LTB spacetimes are the expansion scalar 0, the Ricci scalar 
of the hypersurfaces orthogonal to u", the shear tensor a ab and the electric Weyl tensor E ab : 

~ „ 2R R' , 2\(1-3?)R]' 

o ab = V {a u b) -{®/3)h ab =T.Z ab , E ab = u c u d C abcd = £S ab , (7) 

where h ah = u a u b - g ub , V a = h b V h , and C abcd is the Weyl tensor, Z ab = h ab - 3r]"i] b with 77" = Vh^d? being the 
unit radial vector orthogonal to u". The scalars S and E in (|7| are 



1 

E =3 



R R' 

R~ ~R 



K AT 

^ = -^P + ^. (8) 



The dynamics of LTB spacetimes can be fully characterized by the local covariant scalars {ji, 0, E, $ , 3 &}. Given 
the covariant "1+3" slicing afforded by u", the evolution of the models can be completely determined by a "fluid flow" 
description of scalar evolution equations for these scalars (as in El): 

© = -f-fp-6E 2 , (9) 
P = -P®, (10) 
E = -^E + E 2 -^, (11) 

£ = -|pE-3^f| + Ej, (12) 



together with the spacelike constraints 



(13) 



and the Friedman equation (or "Hamiltonian" constraint) 

2 



K ^ - 



(14) 



The solutions of system 14 1 are equivalent to the solution of the field plus conservation equations Vi } T ab = 0. 



PROPER VOLUME AVERAGE: BUCHERT'S FORMALISM. 

The time slicing defined by u" defines as the space slices the hypersurfaces marked by constant t, whose metric 
is h a b = Uaiib + gui, and their proper volume element is 

dV p = v /det(/i a ,,) dr d0 d</> = .^- l R 2 R' sin 6 dr d0 d</> . (15) 

Consider spherical comoving regions of the form 

&[r] = &[r] xS 2 C 3r [i \ 0[r] = {xeR\0<x<r} (16) 

where S 2 is the unit 2-sphere and x = marks a symmetry center. We introduce now the following definitions: 
Let Z(^[r]) be the set of all smooth integrable scalar functions in 3, then 

Definition 1. For every A £X(@[r]), the p-map is defined as 

where Jq ..dx — J^Zq ..dx. The scalar functions A p : $ — * R that are images of ^ p will be denoted by 
"p-functions". In particular, we will call A p the p-dual of A. 

Definition 2. For every A eX(^), the proper volume average is the functional 

() p [r}:X(®)^R, (A) p [r]=A p (r). (18) 

The real number (A) p [r], associated to the full domain &\f\, will be denoted by the proper volume average 
of A on { 2>\r\. In order to simplify notation, we will drop the "[r]" and p symbols and express ( 18 i simply 



as (A), as it is clear that the average is evaluated as a proper volume integral and it depends on the domain 
boundary as a functional. 

The proper volume average (as well as the p-functions) satisfy the following commutation rule: 

(A)-(A) = (0A)-(0)(A), (19) 

It is important to emphasize that only the functionals (A) can be considered as average distributions, as they satisfy 
(A)) = (A), and thus 

((A- (A)) 2 ) = (A 2 ) -(A) 2 , (20) 
((A- <A)) (!?-<*») = (AB)-(A)(B), (21) 

which define variance and covariance moment definitions for continuous random variables. It is straightforward to 
verify that the functions A p do not satisfy these equations. The relation between the p-functions A p and the average 
(A) is illustrated by figure 1 . 




FIGURE 1. The difference between A p and (A). The figure displays the radial profile of a scalar function A(x) (solid curve) 
along a regular hypersurface £/~^(t), together with its dual p-function A p (x) (dotted curve) defined by jl7| . Panels (a) and (b) 
respectively display the cases when A' < ("clump") and A' > ("void"). The average functional {18} assigns the real number 
(A) [r] to the full domain (shaded area) marked by r) [r] = {x 1 < x < r}, whereas the function A p varies along this domain. Hence, 
A p and (A) are only equal at the domain boundary x = r, and so they satisfy the same differentiation rules locally, i.e. A„(r) = (A) '\r\ 
and A' p {r) = (A)'[r], but behave differently when integrated along the domain. Notice that, from 

then A — (A) < and the opposite situation occurs if A' > 0. This figure also applies for the quasi-local functions and averages. 



46 1 and 447}, if A' < in all & [r] 



The well known evolution equations of Buchert's formalism |4| follow by applying the proper volume average 
functional (18i to both sides of the energy balance (10 1, Raychaudhuri |9]) and Friedman ( fT4] i equations, and then 
using ( 19 1 and (20 1-(21 1 to eliminate averages (A) in terms of derivatives of averages (A)' and squares of averages as 
averages of squares. The averaged forms of these equations are: 

<p+p©> = (p)-+(p>(©)=0, 



(0)2 

3 



= -§<p> 



(&)- 



K 



-3<P> - - 



where the kinematic "back-reaction" term, =2, is given by 



i2[r]^-((©-<©)) 2 )-6<£ 2 ) 



(22) 
(23) 

(24) 
(25) 



Equation (22i simply expresses the compatibility between the averaging ( 18 i and the conservation of rest mass, but 



(23 i and (24i lead to an interesting re-interpretation of the dynamics because of the presence of £2. This follows by 
re-writing these equations as 



(©>• 



(®) 



3R 



cffj 



(26) 



where the "effective" density and pressure are 

ICPeff = K(p)- 



( 3 ^) + £ 



KP e 



cff : 



3 



(28) 



The compatibility condition between (23 1, (24 1 and (25 1 is given by the following relation between J? and ( 3 if )' (as in 
HUCCS): 



J + 2(0) 3 + - (0) ( 3 ^) + ( 3 ^) • 



0. 



(29) 



which is equivalent to the compatibility between the time derivative of ( 14 1, equations (|9|-([T0|, the commutation rule 
( 19 1 and the variance (20 1 for © and From (26 1, the condition for an "effective" cosmic acceleration mimicking 
dark energy is 



J2<0, 



(30) 



2 Lr J 2 

which, apparently, could be possible to fulfill for a sufficiently large and positive back reaction i2. We will evaluate 
this condition for spherically symmetric LTB dust solutions (see [ 16| and lfT2l[T3l for previous work on this). 



QUASI-LOCAL (QL) VARIABLES. 

The Misner-Sharp quasi-local mass-energy function, ./#, is a well known invariant in spherically symmetric space- 
times fT8l[T9l . For LTB dust models ([T]l-(j2]i it satisfies the equations 

UC = KpR 2 R', 2^d = 0, =>Jt = M{r) (31) 

where M is the function appearing in the field equation Comparing ( [31} and |4]) suggest obtaining an integral 
expression for jtf, that can be related to R and R. This integral along the exists and is bounded if we consider an 



integration domain of the form ( 16 1 containing a symmetry center IT91 . Since j$(ct , 0) = for all f, we integrate both 
sides of Q and also ( 3 1 1. This allows us to define a scalar p q as 

K _ 2J( _2M _K \qPR 2 R'Ax 

3 Pl= -j? - -£i - 3 frRlKfo ' (32) 

This integral definition of p q , which is related to p and to the quasi-local mass-energy function, J£, motivates us to 
generalize it to other scalars by means of the following: 

Definition 3: Quasi-local (QL) scalar map. Let X(@) be the set of all smooth integrable scalar functions in 
For every AgX(f), the quasi-local map is defined as 

The scalar functions A q : <2) — > R that are images of J! q will be denoted by "quasi-local" (QL) scalars. In 
particular, we will call A q the QL dual of A. Notice that a quasi-local average can also be defined by means 
of a functional with the correspondence rule (f3"3"]>, but we will not need it in this article. 



Applying the map (33 i to the scalars © and 3 if in (KM) we obtain 



Applying now ( |33] l to p, comparing with ([3]l and Q, and using ( 34 1, these two field equations yield 



® q y k 3 ^ q 



3 J 3- 6' (35) 

% = -^r-2 Pc >- (36) 

Pq = -Pq%- (37) 



which are identical to the Friedman, Raychaudhuri and energy balance equations for dust FLRW cosmologies, but 
given among QL scalars (notice that we are using here the QL functions, not QL averages). By applying ( 33 i to (|8]> 
the remaining covariant scalars {£, #} can be expressed as deviations or fluctuations of p and © with respect to their 
QL duals: 

£=-![©-©,], # = _|[p_p 9 ]. (38) 



SUFFICIENT CONDITIONS FOR A POSITIVE BACK-REACTION AND AN 

EFFECTIVE ACCELERATION. 



Equation ( 30 » provides the relation between back-reaction (Js?) and the an effective acceleration mimicking dark 
matter. A necessary (but not sufficient) condition for the existence of this acceleration in a given comoving domain 
H$[r] of the form ( 16 1 is evidently 



:2={'Hx,r)) >0, 



•fix, r) = (©(*) - (0) [r] f - (&(x) - ® q (x)) 2 , 



(39) 



where we have used (38 to eliminate £ in terms of — & q . Testing the fulfillment of this condition from the integral 
definitions ( 17 i and (33 i is very difficult without resorting to numerical methods. However, for every domain H$[r] 
and every scalar we have A{x) > Vjc S #[r] (A)[r) > (though the converse is not necessarily true). Hence, a 



sufficient condition for the fulfillment of ( 39 1 in a given domain ( 16 1 is simply 

%f(x,r) >0. 



(40) 



Moreover, this condition is still too difficult to handle because of the dependence of on points inside (x < r) and in 
the boundary (x = r) of the domain. Fortunately, by means of the following lemma we can find a condition equivalent 



to ( 40 1 that depends only on the domain boundary and is applicable to any domain (r variable). 
Lemma 1 : (3 s ) = in every domain #[r] for = &(x, r) given by 

P(x, r) = [©(*) - (0) [rf - [®(x) - ® p (x)} 2 , 
Proof . Expanding ( |4T) and applying ( [17) ) and ( [T8| we obtain with the help of ( |20| 



(^)[r] = -(0)[r] : 



%{r) k 



[2@® p -® 2 p }y;dx. 



(41) 



(42) 



Inserting = fp/'fi and © p = i^p/Yp in the integrand above, and bearing in mind that (0) and & p coincide 
at the domain boundary x = r, leads to the desired result: 



(^)[r] = -(0) 2 M- 



1 



T p (r) Jo 



[f p 2 /y p }'dx=-(®) 2 [r}+®l(r)=0. 



(43) 



An analogous result follows for the quasi-local average acting on a scalar like 8? with ( ) q and (/ instead 
of ( ) and 0„. 

As a consequence of this lema, we have ((0(x) — (0)[r]) 2 ) = ((®(r) — ® p (r)) 2 ) and so the sufficient condition for 



J? > given by ( 40 1 can be rewritten now as 



V( r) = ( 0(r ) _ &p{r )f - (0(r) - 0,(r)) 2 > 0, 



(44) 



holding for comoving domains of the form (p~6]>. This condition is domain dependent, in the sense that it may hold for 



some domains and not for others. Since the condition for an effective acceleration in ( 30 1 is given as the average of the 
scalar (2/3)^ '(x, r) — (k/2)p(x), lemma 1 also provides the following sufficient condition for its fulfillment 



^ {r )= cf (r) _ -p(r)>0. 



(45) 



where c (o(r) is given by (44 1. For the remaining of this paper we will examine conditions (44i and (45 i, finding out 



first the necessary restrictions on LTB models to comply with (|44J, and then exploring for these models the fulfillment 
of d45l. 



PROBING THE SIGN OF THE BACK-REACTION TERM. 



In order to look at conditions ( |44) > and ( |45] l, we need to explore the behavior of LTB scalars along radial rays of 
the hypersurfaces £?( 3 K For this purpose, the integral definitions (17i and (33 1 yield the following properties of p- 
functions and averages and their quasi-local analogues: 



, 3R' & 



A'^^iA-A,} 



A(r)-A p (r) = ^ A>(x)Y p (x)dx, A(r) - A q (r) = jf A' '(*) f q (x) dx, 



(46) 
(47) 



where & p is the p-function associated to &(r) and ^ = / & dY p = An^R 2 R'dx = (47t/3)R 3 (r). Considering (47 1 



condition d44l for £2 > becomes 



(r) = #(r)¥(r) > 0, 



and with <p and y/ given by 



with: 4>(7) = / &'(x) <p(x,r) dx, 
Jo 



¥(r)= / 0'(x)y/(x,r)dx, 
Jo 



(p(x,r) = 
y/(>,r) 



W W y p (r) 



^p(r) 



(48) 

(49) 
(50) 



where we have used the relation 'Vqj'Vp — & p , which follows directly from ( 17 1 and (33 i and is valid for all domains. 

The fulfillment of ( 48 1 clearly depends on the signs of <p and y/ (besides the sign of ©') at all points in any 
given domain. This fulfillment might hold only in some domains and not in others. Since, by their definition, 
f p (0) = ig(0) = and #"(0) = 1, ( 50 1 and d49l imply that for every domain we have 1/(0, r) = 0, y/(r, r) = 2, whereas 
<p(0,r) = <p(r : r) = 0. Thus, as long as & and/?' are non-negative for all jc £ #[r], the sign of y/ is non-negative, and 
so *P basically depends only on the sign of ©'. On the other hand, the sign of <p is not determined, and so the sign of <t> 
requires more examination as it depends on both: the sign of ©' and the ratio 3? p {x) / ' & p (r) (which relates to the sign 
of Since & and R' can become negative in elliptic configurations where the ,3^^\t) have spherical topology, we 
will only consider in this paper "open" LTB models in which these slices are homeomorphic to R 3 . Since the sign of 
'tf depends on the sign of the product O^P, we will need to obtain the conditions for both terms having the same sign. 



In order to to probe condition (48 i, we consider the following sign relations that emerge from (46 1— ( 47 i and are 
valid for all x S # [r] in every domain and in every 

^(3). 



jT'<0 



J?(r)>& p (x), 



(51) 
(52) 



Considering (48 i— ( 50 1 and the relations above, we have the following rigorous results: 



Lemma 2 : Let © be a monotonous function in a domain ( 16 1 in a given 

^(3) 

, then 



4[r] = 
£\r\ > 
2[r] < 



if the domain belongs to a parabolic model, 

if the domain belongs to a hyperbolic model, 

if the domain belongs to an open elliptic model for which 



<0, 



(53) 
(54) 
(55) 



The proof follows directly from (48 1— ( 50 1. It is trivial for the parabolic case, since J? = 1 for every domain. 



For the hyperbolic and elliptic cases, having assumed a monotonous © implies that the sign of c £ only 
depends on the sign of (p(x,r) in (49 1, but then regularity conditions ITOlfTTI require > to hold for all 



r in hyperbolic models, while regularity admits (but does does not require) that < holds for all r in 
open elliptic models (we are not considering closed models, see [ 16 1 for the study of back-reaction in these 
models). Hence, we have > and < 0, respectively, for the hyperbolic and elliptic cases, and the result 
follows. 



Notice that (53 i is domain independent (holds at every domain #[r]), whereas (54i necessarily becomes domain 



independent in those of hyperbolic models for which © is monotonous for all r. For hypersurfaces of 



open elliptic models in which © is monotonous for all r, then ( 54 1 may become domain dependent: it holds for all 
domains sufficiently close to the center (where < 0) but may not hold if changes sign, though ( 54 1 is domain 
independent if does not change sign (see lfl6l ). 

The restriction that © be monotonous for all r does not hold, in general, for all in any hyperbolic or open 
elliptic model (see [20] for details). However, the effect of 0' changing sign simply makes the sign of Js? domain 
dependent along any given i^( 3 ). The same remark holds for open elliptic models for which changes sign: 

Lemma 3 : Consider the following two possible configurations along a given 5^ 3 ': (i) ©' chan ges sign at a 
given i = &[r] in a hyperbolic model, and (ii) and 0' respectively change sign at a x — y\ , x — yi 
inside #[r] in an open elliptic model. For both cases there are always domains #[rn] such that Js?[r] > 
holds for all r > tq. 

The proof in both cases is based on the fact that the zero of ©' is fixed at any given SF^\ whereas the choice 
of domain is arbitrary. Consider the hyperbolic case: the zero of 0' implies that © is monotonous in the full 
"external" range r > x > y. Hence if regularity conditions hold, r can take arbitrarily large values and we 



can always find sufficiently large domains such that r ^> y holds, and thus we can use ( 5 1 1 to obtain the 



expected sign in the integration of <t> and VP in (48 1 along the range y < x < r. For sufficiently large domains 
the "external" contribution y < x < r outweighs that of the inner range < x < y. The same situation occurs 
in the elliptic case, for which must pass from negative to positive at x — y\. The zeroes of 0' and 
imply that F and are monotonous forx > max(ji,y2), an d thus the same argument applies. The reader is 
advised to consult [ 16] for further detail on these proofs. 

As a consequence of lemmas 2 and 3, the condition £2 > is compatible with both, negative and positive spatial 
curvature (hyperbolic and open elliptic models), at least in their asymptotic radial range. We need now to find if the 
positive back-reaction is comparable to the average density. 




FIGURE 2. Scenarios with positive effective acceleration. The figure displays, for a density clump (a) and a void (b), the radial 
profiles of p(x) (black solid curve) and ^(x) (red solid curve), compared with (p) (black dotted curve) and back-reaction J2 oe (^) 
(red dotted curve) along a given regular hypersurface ST^' of either a hyperbolic or an open elliptic model. Notice an intermediate 
region (gray area) where density is low and nearly homogeneous, while the growth of back-reaction is driven by the gradients of 
spatial curvature and reaches a maximum, thus allowing for |45j to hold. 



PROBING THE SIGN OF THE EFFECTIVE ACCELERATION. 



In order to examine ( |45) , we will need the following rigorous results concerning the radial profiles of scalars A along 
the 

Lemma 4 . If R' > everywhere and A (7) — ► Ao = constant as r — ► °°, then A p (r) — > Ao and A 3 (r) — > Ao as 
r — ► oo, while A p (r) — > and A^(r) — > in this limit. 

Proof . We consider the case when A' > for sufficiently large values of r. The case when A' < is 
analogous. If the limit of A as r — > °° is Ao, then for all £ > there exists L(e) such that Ao — £ < A(r) < Ao 
holds for all r > L(e). Constraining A by means of this inequality in the definitions of A p and A q in ( 17 1 and 
(33 i leads immediately to Ao — £ < A q (r) < Aq and Ao 
r — > oo is also Aq. The limits A' p (r) — » and A' (r) 



£ < Ap(r) < Ao. Hence, the limit of A p and A q as 
follow trivially. 



Lemma 5 . If there is a zero of A' at x ■ 
of A p at x — r\ > y and a zero of A^ at x 



y and /?' > in #[r], then for sufficiently large r there will be a zero 
r 2 >y, withA(ri) =A p (r 1 ) and A(r 2 ) =A q (r 2 ). 



Proof. Let A' pass from positive to negative at x — y. As x reaches y the first integral in ( |47| > is still positive 
and so A(y) >A p (y), but for y < x < r the integrand becomes negative, and so the contributions to the integral 
are increasingly negative. Since "V p {x) /^p(r) is increasing, if r is sufficiently large, then a value x = r\ > y 
is necessarily reached so that this integral vanishes (thus A(n) = A p (n)). From (46 1, we have A' (ri) = 
and A p < for* > r\. An analogous situation occurs when A' passes from negative to positive. The proof is 
identical for A q , but using the second integral in (47 1. 



In order to apply these lemmas to the case A 
the gradients of & q and @ p as 



with c & given by (44 1, we use (46 1 to rewrite this scalar in terms of 



R 

3R' 



@ 



/ 2 



R 

3R> 



©', 



J_ 



0' 



(56) 



Considering now LTB models (hyperbolic and open elliptic) in which the scalars {0,p, 3 ^} tend to nonzero finite 
values {©o,po, 3 ^o} as r — ► °° ( an asymptotic FLRW state in the radial direction), then Lemma 4 and (56i imply that 
^€ — * 0, — ©„ — * and — O — ► as r — * °°, and also 0'^ and 0^ vanish in this limit. Since tf(0) = follows 



from (47 1, then domains must exist for which *€' must have a zero for some x = y £ #[r]. Lemma 5 implies then that 
domains must exist in which £1' = (2/3) ' has also a zero at ri > y, corresponding to a local maximum of £2 where 
it reaches its maximal value in the domain (see figure 2). 



The next step is to compare c i and i? with p and (p) in order to test the fulfillment of (45 I. For this purpose we note 
that 0, p , are respectively related to p, p p , p q and 3 ^ p , ^S% q by the constraints ( 14 1, (24 1 and (35 1. Hence, 
the magnitude of the expansion gradients is closely connected with the magnitude of the gradients of the density and 
spatial curvature. Given the fact that =S(0) = and i? — > as r — > °°, reaching a maximum in the intermediate range, 



the best possible situation in which pB} could hold is if £2 reaches its maximum in the same region where p (and 
thus (p)) has a low and almost constant value. In this scenario we would have a large intermediate region in which 
p(7)/p(0) < 1 and p' sa (and thus (p)[r]/p(0) < 1 and (p)' w hold), so that the growth of c i (and thus J2) is 
driven by the gradients of the spatial curvature. The scenario described above is illustrated in figure 2 for a density 
clump and void profile. 



CONCLUSION. 

We have introduced quasi-local variables and averages in LTB dust models. These scalar variables are very useful 
to discuss various theoretical issues concerning these models, such as the application of Buchert's scalar averaging 
formalism JH Q. We have found in this paper analytic conditions for a positive back-reaction term £? and for 
an effective acceleration mimicking dark energy in these models. The present paper provides a quick summary of 
comprehensive articles |[T5l [T6 20 1 that are currently under revision, all of them dealing with different aspects of radial 
profiles of scalars and the issue of back-reaction in LTB models. We have chosen here the simplest boundary conditions 



(radial asymptotical homogeneity) to illustrate how generic LTB configurations with open topology (hyperbolic and 
elliptic) can fulfill the conditions for such an effective acceleration. More general conditions are examined in 11201 . It 
is very likely that the astrophysical and cosmological effects of Buchert's formalism will require numerical methods 
applied to more "realistic" configurations not restricted by spherical symmetry and compatible with observations. 
The analytic study carried on here and in the associated references can certainly provide a useful guideline for this 
important task. 
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